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LINEAR  FILTERING  FOR  SHIPBORNE  MAD  RECORDERS 


Summary 


This  report  describes  the  derivation  of  linear  filters  for  MAD 
equipment  to  be  used  for  submarine  detection  from  high-speed  surface 
vessels.  Since  the  bandwidth  required  to  pass  all  the  anomaly  signals 
of  interest  is  considerably  larger  than  in  airborne  systems,  it  is  recom- 
mended that  the  passband  be  covered  by  two  filters,  one  driving  each 
pen  of  a two-channel  recorder.  It  is  also  recommended  that  each  of 
these  filters  be  provided  with  switchable  cutoff  frequencies  for  improved 
noise  discrimination  whenever  the  target's  closing  rate  can  be  estimated. 
Furthermore,  a selection  of  variable-frequency  notch  filters  should  be 
available  to  cascade  with  the  basic  filters  in  order  to  remove  quasi- 
periodic  noise,  both  that  of  geomagnetic  origin  as  well  as  periodic 


maneuver  noise. 
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1.  Introduction 

The  class  of  signals  to  bo  considered  is  that  gen- 
erated by  passing  through  the  field  of  a magnetic  dipole  with 
speeds-along-track,  v,  between  20  and  70  knots  and  smallest 
slant  ranges,  z,  between  300  and  1000  feet  (figure  1). 


Additive  noise  corrupts  the  observed  signal.  This  noise  is 
assumed  to  be  the  sum  of  a stationary  Gaussian  component  with 
power  spectral  density 

N(co)  « N Q/J-  (1.1) 

within  the  passband  of  the  signal  collection,  and  a nonstationary 
component  comprising  a random  sequence  of  quasi-periodic  wave- 
trains  caused  both  by  natural  geomagnetic  phenomena  (TH,  pp. 

7-55#  [2],  p.  287£f0  and  by  periodic  maneuvers  (trapping  circles, 
cloverleaves,  etc.  ([ll#  PP.  67-71)). 

Filters  are  required  to  discriminate  against  the 
noise  while  preserving  the  shape  (in  time)  of  the  anomaly  signals 
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as  closely  as  possible.  Although  the  center  frequency  of 
naturally  occurring  quasl-periodic  noise  Is  relatively  stable 
once  the  wavetraln  has  begun,  this  frequency  Is  unpredictable 
with  any  precision  In  advance  of  the  wavetraln' s commencement; 
similarly,  the  center  frequency  of  periodic  maneuver  noise  will 
presumably  be  determined  by  as  yet  unknown  tactical  consider- 
ations. As  these  quasi-periodic  phenomena  have  extremely 
narrow  bandwidths  with  respect  to  the  bandwidth  of  any  anomaly 
signal,  it  is  evident  that  a nearly  optimal  linear  way  of 
minimizing  their  effect  is  with  an  adjustable  notch  filter 
tuned  with  its  frequency  of  infinite  rejection  at  the 

center  freouency  of  the  disturbance.  It  is  therefore  suggested 
that  at  least  one  such  tunable  notch  filter  be  available  to 
prefilter  the  signals  before  further  processing. 

The  remainder  of  this  report  is  concerned  with  fil- 
tering against  the  stationary  Gaussian  component  of  the  noise. 
Since  the  bandwidth  of  this  noise  is  large  compared  to  the 
bandwidth  of  an  anomaly  signal,  it  oannot  simply  be  stripped 
away  with  a narrow-band  filter.  It  is  intuitively  clear  that 
the  optimal  linear  filter  will  be  one  whose  response  at  any 
frequency  is  directly  related  to  the  ratio  of  signal  power  to 
noise  power  at  that  frequency;  this  notion  will  be  made  more 
precise  later  on.  It  follows  immediately,  however,  that  sinoe 
the  noise  power  is  not  equally  distributed  among  all  freauencies 
(equation  (1.1))  it  is  unlikely  that  the  filter  will  have  con- 
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stant  gain  over  the  passband  of  an  anomaly  signal.  Therefore 
some  waveform  degradation  is  inevitable;  indeed,  it  can  be 
shown  that  the  optimal  linear  filter  for  a single  signal  in 
noise  with  spectral  density  given  by  equation  (1.1)  is  a 
"matched  filter"  for  the  signal's  derivative,  whose  output(in 
the  absence  of  noise)  is  the  first  derivative  of  t he  auto- 
correlation function  of  the  signal.  Conversely,  a filter  that 
preserves  signal  waveshape  exactly  can  remove  no  noise  at  all. 

In  the  sections  of  the  report  to  follow,  filters 
will  be  designed  based  on  the  theory  of  optimal  linear  fil- 
tering due  to  N.  Wiener  and  A.  N.  Kolmogorov.  It  will 
unfortunately  turn  out  that  the  filters  so  specified  are 
unrealizable  with  the  customary  physical  components,  and  must 
be  approximated.  In  the  interests  of  simplicity,  ease  of 
actual  synthesis,  and  freedom  from  critical  adjustment,  it  was 
decided  to  make  the  approximation  with  physically  realizable 
filters  having  only  real  poles  and  zeros. 

2.  Development  of  a mathematical  model 

It  has  been  shown  (113,  p.  1^7  ) that  if  the  vector 
closing  velocity  v (figure  1)  is  constant,  than  any  anomaly 
signal  can  be  represented  as  a linear  combination  of  three 
linearly  Independent  basis  signals 

s(t)  * a0s0(t|  + a1s1(t)  + a2s2(t) 
in  which  the  coefficients  deP0nd  ln  a complicated  fashion 
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on  geographical  location,  dipole  orientation  and  strength,  and 
other  time -independent  physical  parameters.  The  three  basis 
signals  are  given  by 

e0(t)  = (02+l)‘5/2 

8l(t)  = 

s2(t)  = (e2+i)‘5/2e2  , 
where  0 is  a dimensionless  variable 


0 = t/c 

and  t is  a parameter  with  the  dimensions  of  time 

X «=  z/v  = (smallest  slant  range  )/(c  losing  speed). 

In  section  1,  it  was  given  that  300  <:  s ^ lf?00  ft.  and 

20  £ v £ 70  kn.;  for  the  surface  vessel  MAD  problem,  then,  the 

parameter  x:  must  lie  between 


and 


rmin  * = 

v 

max 


300  ft  = 2.5U  sec. 

70kn.  x 1 ,687^f t/sec/kn 


max 


'max 


1500  ft 


vmln  20kn.  x 1 .6fl7Bft/secAn 


= U4.ll  seo. 


(2.2) 

(2.3) 

(2.U) 


(2.5) 

(2.6) 


For  interpretation,  these  figures  mean  that  the  bell-shaped 

signal  $Q(t),  equation  (2.2),  will  occur  with  pulse  widths 

between  2.86  sec.  (when  Z B f-*.,)  and  5>0.2  sec.  (when  X « Z ) 

mj.n  max 

at  half  amplitude,  or  between  3.6U  and  63<6  sec.  at  10#  amp- 
litude. For  a second  interpretive  example,  "W"  signals,  s2(t), 
will  be  observed  with  a "time  between  peaks"  of  from  U.15 
( C=Tmin)  fco  72 *5  sec.  (^=Cmax). 

Thus  an  anomaly  signal  is  completely  described  by 
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stating  5 numbers : aQ,  a^  and  a?  which  govern  the  signal's 
shape  by  equation  (2.1),  Z,  which  sets  the  time  scale  of  the 
signal,  and  t^,  the  "epoch  instant"  or  time  of  occurrence, 
corresponding  to  t=0  in  equations  (2.2).  But  since  the  purpose 
of  MAD  is  the  detection  and  localization  of  suspected  anomalies, 
it  is  unlikely  that  any  of  these  5 quantities  can  be  known  in 
advance.  It  is  not  unreasonable,  therefore,  to  describe  them 
as  random  variables  with  an  unknown,  and  perhaps  time -varying. 
Joint  probability  distribution.  Furthermore,  in  a tactical 
encounter  it  is  to  be  expected  that  a sequence  of  anomaly 
signals  with  differing  values  of  ^aj^,  Z,  and  t^  will  be 
observed  occurring  in  bursts  or  clumps,  corresponding  to 
repeated  detections  of  the  ?ame  anomaly. 

If  the  five  parameters  describing  a single  signal  are 
regarded  as  random  varaibles,  then  a sequence  of  signals  may  be 
treated  as  a random  process.  This  process  is  likely  to  be  non- 
stationary, for  consider  a successful  encounter  culminating  in 
the  licalization  of  a stationary  faarget:  as  the  encounter  pro- 
ceeds in  time,  the  quantity  T will  (on  the  average)  decrease 
on  successive  detections  as  the  range,  z,  decreases.  Including 
nonstationarity  in  the  description  of  the  signal  process, 
however,  must  ultimately  result  in  the  specification  of  a time- 
varying  filter.  For  the  sake  of  simplicity  in  computation  and 
eventual  filter  construction,  the  signal  process  will  be  taken 
to  be  stationary;  it  is  assumed  that  nonstationarity  (like  that 
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in  the  mythical  encounter  described  above)  can  be  adequately 
treated  in  practice  by  making  some  filter  parameters  switch- 
able  by  the  MAD  equipment  operator.  Sufficient  conditions  for 
stationarity  are  (1)  the  quantities  an<*  Z ^ave  a 

Joint  distribution  and  their  values  on  successive  anomaly 
signals  are  statistically  independent,  and  (2)  the  epochs  t^ 
are  uniformly  and  independently  distributed  over  any  finite 
interva 1 . 


For  the  filter  design  theory  in  the  next  section,  the 
signal  process  is  sufficiently  known  when  its  x$h  power  spectral 
density  is  given.  Since  the  epochs  are  uniform  and  independent, 
Campbell’s  theorem  ( C 33  » P-  ) is  applicable,  with  the 

result  that  the  power  spectral  density  of  the  process  is  pro- 
portional, at  each  frequency,  to  the  energy  spectral  density 
of  a single  anomaly  signal.  Taking  the  constant  of  proportion* 
ality  equal  to  1 for  convenience,  and  letting  F be  the  joint 
c.d.f.  (cumulative  distribution  function)  of  the  random  variables 
|aj|p  and  Z,  the  power  spectral  density  S(oo)  of  the  signal 
process  is  given  by 

S(oj)  = JdF(a0,  a1#  a?,  T)  |7{s(t)|}2  (2.7) 

where  ^ indicates  the  operation  of  Fourier  transformation, 
and  s(fc)  is  given  by  equation  (2.1).  If  the  additional 
assumption  is  made  that  the  aj  are  pair-wise  uncorrelated,  the  / 

cross-terms  in  the  integrand  of  equation  (2.7)  disappear  and 
S(io)  can  be  written 
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S(*0  = 

where  G is  the  c.d.f. 

3.  Filter  design  - theory 

With  the  signal  and  noise  models  thus  far  developed, 
the  data -processing  problem  clearly  belongs  in  the  domain  of 
the  Wiener-Kolmogorov  theory  of  optimal  smoothing.  An  obser- 
vation is  made  of  the  stationary  random  process  v(t)  with 
spectral  density  (from  eauations  (1.1)  and  (2.8)) 

V (<o)  = S(to)  + Nq/<o2 

(assuming  signal  and  noise  are  at  least  uncorrelated),  and  as 
faithful  a copy  of  s(t)  as  possible  is  to  be  extracted.  The 
minimum-mean-square-error  linear  realizable  zero-lag  smoothing 
filter  H ( Jto)  is  in  principle  found  by  solving  an  appropriate 
Wiener-Kopf  integral  equation;  however,  by  exploiting  the 
relatively  simple  noise  structure  H(j^)  can  be  found  by  simpler 
techniaues.  If  H(Jio)  is  the  best  linear  filter  for  s in  s+n, 
then  (by  linearity)  it  is  also  the  best  linear  filter  for 
s in  s+n.  But  here  the  power  spectral  density  of  n is 
|jujp(NQ^j  ) = Nq;  i.e.,  n is  white  noise,  and  the  best  linear 
filter  for  any  signal  in  white  noise  can  be  written  down  by 
inspection  (Llj.],  p.  72  ff.).  For  the  problem  at  hand,  H(j«) 
is  the  minimum-phase  filter  whose  real  part  is  «*>2s(uj): 

H(ju)  =to2s(u>)  + 

where  denotes  the  operation  of  Hilbert  transformation. 


I a,2  f |s<(w;C)|2  dG(r), 

I r0  J 

of  T,  and  Sj(u);x)  ■ '7{sj(t»'cd* 
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It  is  therefore  necessary  to  return  to  equation  (2.8) 
and  compute  u^S(uO;  the  first  step  is  clearly  to  find  the 
Fourier  transforms  (to;  r)|  2 of  the  signals  3q,  s^  and  Sg. 

Use  of  the  Bateman  tables  (C53*  p.  9 ) yields 

<o2  | SQ  (to;  X)  I 2 = *6K22  ( ru>) 
a;2  |SX (*o;xr)|  2 = a^K12(ru;) 

^2|S2(^;r)|2  = [7u2k1  (rco)  - z^K^r^f]2 
(up  to  a common  factor).  The  Kj(»)  are  modified  Bessel  functions 
of  the  second  kind.  Now  the  right-hatad  sides  of  equations 
(3*3)  should  be  inserted  into  equation  (2,8)  to  compute  ^S(u>). 
Unfortunately,  neither  the  mean-square  coefficient  values 
aj^  q nor  the  time-scale  c.d.f.  G(T)  are  known. 

In  order  to  proceed  with  the  design,  an  ad  hoc  tech- 
nique is  employed  which,  while  by  no  means  optimal,  is  nonetheless 
plausible  and  does  lead  to  a comparitive ly  simple  design.  In 
brief,  io2S(to)  is  taken  to  be  a bandpass  spectrum  with  its  low- 
frequency  behaviour  determined  by  that  of  u?  |sQ(««J;'CmHy  )|  2 and 
its  high-frequency  behaviour  determined  by  that  of 
u>2{so(u);Train)|  2.  (A  rationale  for  this  technique  is  given  in 
Appendix  1 to  this  report^)  With  this  specification,  the  desired 
filter  is  described,  using  equation  (3*2)  by 

fEW rmax ) ^K2  ^ rmax 0 2 t*\  s ..A«-2 

Re|H(J<o)|  * | constant 


(3.3)  j 


2(^rmin^2 


Tmax 


UJ  > 


tmin 

-LlS. 


(34) 


fmax 

The  number  1.9  is  determined  by  the  fact  that  the  function 
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x^K2(x)  assumes  a maximum,  with  slops  zero,  at  x=1.9  (L6]f  p.4  17) 

Of  course,  a finite  RLC  circuit  cannot  be  built  to 

satisfy  equation  (3.4),  since  the  real  part  of  any  RLC  trans- 

2 

mittance  must  be  a rational  function  of  . Therefore  the 
problem  becomes  one  of  approximating  equation  (3*4)  with  the 
real  part  of  an  RLC  transmittance. 


4.  Filter  design  - approximation 

Let  H(s)  be  the  RLC  transmittance  whose  real  part 
approximates  equation  (3.4);  i.e.,  Re |h(  Re|H(Jo>)^. 

Then  H(s)  is  a rational  function  of  s.  Since  it  is  clear 
(from  eauation  (3.4) ) that  Re^H(j«)|  is  never  negative,  it 
may  tentatively  be  assumed  that  the  degrees  of  the  numerator 
and  denominator  of  ft(s)  differ  by  exactly  one.  Thinking  of 
H(s)  as  a high-pass  and  a low-pass  filter  in  cascade,  it  is 
plausible  that  the  high-pass  and  low-pass  sections  may  (to  a 
rough  approximation)  be  designed  independently  of  one  another 
since  1«9/Tmax  and  1.9/rmln  (see  equation  (3*4))  differ  by 
about  1 1/4  decades. 

On  tho  low  end,  lim  x^Kp^(x)  = 4x^  ( L61 , so 

x-*-0 

that  the  minimal  high-pass  section  of  H(s)  is  of  the  form 


As  beoause  lim  Re[  AJw  ? = A . Hence  the  form  of 

+ p0  *»“p,0  + Po»  W 

ppro: 

H ( s ) 


a + p0  ~-*0  PqI  PO 

the  approximating  filter  is,  as  a first  guess. 


As  s + z-i  S + Zo  s + Z3  s + Zl|  1 

s + P0  s + px  s + p2  s + p3  s + p^  s + P£ 


(4.D 
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with  tha  zeros,  poles  and  the  factor  A to  be  chosen  so  as  to 
achieve  the  best  approximation  of  Re^H(j“->)|  to  Re^H(  j«)j  . 

One  way  of  instrumenting  the  approximation  procedure 

A 

is  shown  in  figure  2.  H(s)  is  built  with  all  frequencies 


Re{$(J«)} 


Figure  2 

scaled  up  into  the  high  audio  range  and  with  the  ten  free  para- 
meters of  equation  (4.1)  available  for  adjustment.  The  system 
output  is  displayed  on  an  oscilloscope  against  a transparency 
of  Re{H(juJ)|  and  the  parameters  A,  ^ZjJ^  and  {pi£o  adJu3ted  by 
hand  to  give  a "good"  approximation.  Not  only  are  the  equip- 
ment requirements  to  implement  this  scheme  quite  stringent 
(e.g.,  the  low-pass  filter  must  be  tunable  to  remove  the 
second  harmonic  of  the  oscillator  as  it  is  swept  over  several 
decades),  but  the  adjustments  of  the  ten  parameters  would  be 
interdependent  rendering  a good  approximation  difficult  to  find. 
The  approximation  method  finally  chosen  employed 
digital  computation?  Re|H(Jus)j  was  sampled  at  94  points  over 
the  frequency  range  0 to  approximately  2 c/a,  and  a digital 
computer  was  programmed  to  employ  a modified  Newton's  technique 

A 

that  varied  the  parameters  in  H(s),  starting  from  a set  of 
initial  estimates,  so  as  to  minimize  the  squared  error 

e2  = *|  [Re{H(j«k)}  - Re$H(  j^)}]2. 
lc— 1 


#A  description  of  the  method  is  given  in  Appendix  2. 


(4.2) 
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Several  different  sets  of  initial  parameter  estimates  were 
used,  in  order  to  make  more  likely  the  achievement  of  a global, 
father  than  just  a local,  minimum  for  £ . The  Li.-zero,  5-pole 
filter  that  resulted  was  deemed  to  have  an  unacceptably  large 
squared  error;  furthermore,  its  3 db  bandwidth  was  approx- 
imately 2.2  decades.  Since  the  3 db  bandwidth  of  a single 
bell-shaped  signal  SqH)  is  only  about  0.6  decades,  it  was 
felt  that  a 2.2  decade  wide  filter  would  pass  so  much  noise 
that  weak-signal  performance  would  be  very  poor. 

Sinoe  existing  MAD  equipment  is  provided  with  two 

display  channels,  it  was  decided  to  cover  the  frequency  range 

of  interest  with  two  filters,  the  first  for  signals  having 

values  of  the  time-scale  parameter  x.  between  rrain>HP  = 2. 54 

and  rmax  = 10.6  seconds,  and  the  second  for  the  range 

between  r . Tta  * 10.6  and  t __  = Uii.U  seconds.  Two 
ram  i uv  max  | LF 

filters  having  the  form  of  equation  (4.1)  were  designed  using 
the  computer  program  mentioned  above.  After  some  experimentation. 


it  was  observed  that  a considerable  reduction  in  minimum 


squared  error  was  achieved  by  deleting  the  numerator  factor 
(s  + z^)  from  each  filter.  This  is  reasonable,  in  retrospect, 
since  the  real  part  of  each  filter  function  becomes  negative 
at  the  high  end  of  its  passband  allowing  a better  fit  to 
Re|H(j&>)]  which,  though  always  positive,  decreases  exponentially 
to  zero  at  high  frequencies.  It  was  also  noted  that  the 
minimum  squared  errors  of  the  best  5-pole,  3~zero  filters  were 


1 


CONFIDENTIAL 


11 


ouivtmuitM- 


not  significantly  Increased  by  deleting  a factor  ( s+z^ )/( s+pcj) 
from  each  filter  and  optimizing. 


5.  Filter  construction  and  estimated  performance 

The  designs  finally  arrived  at  are  given  by 

Hlp(8)  ( s + . 0213 ) (I+. 459 )( s+1^92 Hs+:7S7r 

HhpU)  = ( sT;o59Birsll .961  y ( ) ( s+3 .2651  ‘ 

Because  of  the  extremely  low  frequencies  involved,  the  use  of 
inductance  in  the  construction  of  the  filters  is  impractical. 
However,  it  is  well  known  that  any  RLC  transmittance  can  be 
realized  with  RC  circuits  only  in  the  forward  and  feedback 
paths  of  an  ideal  operational  amplifier;  this  method  of  synthesis 
is  certainly  appropriate  here.  The  squared  error  between 
Re[fij  and  RejH]  Is  at  a minimum  when  the  filter  poles  and 
zeros  are  those  given  in  equations  (£.1).  From  a practical 
standpoint,  though,  the  squared  error  is  moderately  insensitive 
to  small  errors  in  the  poles  and  zeros;  *>%  tolerance  is  prob- 
ably adequate.  The  squared  error  is  least  insensitive  to  the 
lowest -frequency  pole  in  the  two  filters  (s  = -.0213,  a * -.0898), 
because  the  noise  power  density  is  highest  at  low  frequencies. 

But  this  fact  is  relatively  unimportant,  since  it  is  recommended 
in  section  6 that  these  poles  be  adjustable  anyway.  Values  of 
the  gain  factors  (1.75  and  30.21;.)  are  of  course  unimportant; 
indeed  it  can  be  shown  that  with  these  gain  factors  and  an 


(5.1) 
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input  noise  spectral  density  Vq/*>  , the  r.m.d.  output  noise 
from  the  HP  filter  is  slightly  less  than  one-half  the  r.m.s. 
output  noise  from  the  LP  filter.  Thfcs  these  gain  factors 
should  be  individually  adjustable.  In  practice , the  gains 
would  probably  be  set  so  that  the  r.m.s.  noise  outputs  from 
the  two  channels  were  approximately  equal. 

An  alternative  synthesis  is  Implied  by  expanding 
equations  (£.1)  in  partial  fractions: 


Hlj?(s)  - -.107?  + 4.9378  + i,2g9< 
s+,0213  s+,459  S+.59J 


8m(.) 


s+7099 


s+.787 


+ -28.1646  . 
s+3 .206 


Realization  according  to  equations  (£.8)  would  require  two 
operational  amplifiers  for  each  filter;  each  Operational 
amplifier  would  have  two  simple  RC  Integrator  networks  con- 
nected to  its  summing  junction. 

Figure  3 i*  « graph  of  RejHj  and  Re^H^  for  the  two 
filters.  In  each  case  the  squared  error  is  £2  *»  8.1.  Since 
the  sum  of  the  squares  of  the  94  sample  points  of  Re[H^  is  833* 
the  percent  r.m.s.  error  per  sample  point  is 

A **  ’ 45?  x 100  * 10<‘ 

That  is.  Re^H^pl  and  match  the  desired  response 

(equation  (3*4))  within  an  r.m.s.  error  of  lo£  at  eaoh  sample 
point . 

Figure  4 shows  Jhlp(J«*»)|  and  |SHP(j~)| , the  gain 
functions  of  the  two  filters.  Their  (oommon)  3 db  bandwidth 
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(5.3) 
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la  approximately  1.3  decades. 

Figures  5 and  6 are  computed  waveforms  indicative  of 
the  general  form  of  the  output  signals  to  be  expected  from  the 
filters.  (Appendix  3 gives  a typical  filtering  algorithm  used.) 
Figure  £ shows  the  output  from  the  LF  filter  when  driven  by 
the  slowest  signal  designed  for:  sQ(t,  44.4) • Since  this  input 
signal  falls  at  the  low  edge  of  the  passband,  the  filter  does 
some  differentiation  and  causes  the  peak  of  the  output  signal 
to  occur  before  the  peak  of  the  input  signal.  Figure  6 shows 
the  output  of  the  LF  filter  with  an  input  of  Sg(t;  10.6),  the 
fastest  signal  it  is  designed  to  handle.  In  this  oase  the 
filter  does  some  integration,  and  the  output  peak  is  delayed 
with  respect  to  the  Input  peak. 

Figures  7 and  8 extract  two  of  the  most  informative 
features  of  the  filters'  output  waveshapes,  and  show  them 
plotted  against  the  input  signal  time -scale  parameter,  t . The 
input  signal  for  these  figures  is  the  bell-shaped  signal  SgftjT), 
which  (for  every  t)  has  peak  value  1 at  t=0.  It  was  noted 
above  that  the  time  between  peak  input  signal  and  peak  output 
signal  depends  on  the  time  scale  of  the  inptt  signal;  the 
curves  labeled  "TIME”  display  this  dependence  for  both  filters. 

In  view  of  the  dependence,  it  is  clear  that  some  care  must  be 
exercised  if  detected  anomaly  signals  are  to  be  used  to  estimate 
the  direction  of  an  anomaly. 

The  curves  labeled  "PEAK"  in  figures  7 and  8 can  be 
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Fig.  7 HF  Filter  - Peak  Output  for  IV  Peak  Input 
and  Time  from  Input  Peak  to  Output  Peak 
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Fig.  8 LF  Filter  - Peak  Output  for  lV  Peak  Input 
and  Time  from  Input  Peak  to  Output  Peak 
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uaod  to  estimate  the  overall  performance  of  the  filter  In  a 
noisy  environment.  Suppose,  as  a reasonable  criterion,  that 
an  anomaly  signal  is  detectable  If  at  the  output  of  the  filter 
its  peak  value  is  equal  to  the  r.m.s.  noise  level.  Then,  for 
example,  if  the  r.m.s.  noise  level  at  the  output  of  the  LP 
filter  is  1 If,  an  anomaly  signal  Sq( t jUU .l|. ) of  peak  amplitude 
(see  figure  0 at  TT=44.4)  1/2.02  = .354V  is  detectable. 
Pursuing  the  example  further,  the  r.m.s.  output  noise  is 

nrms  * (_/nm|Hlj,(J~)|2  d^TT  )*, 
where  N(«o)  = NQ/<c£,  equation  (1.1),  and  HLp(s)  is  given  in 
equation  (5  1 ) • The  value  of  the  expression 

* 

( J (1/*52)|hlp(J^)|2  d*3/2ir  )® 

• <70 

is  readily  computed  to  be  22.15,  so  that 

nrms  = 22.15 /N^. 

Hence  if  n = 1 V then  N = (1/22/15)2  = .00204;  at  a fre- 
rms  o 

quency  of  .01  c/s  (<*>=.1256),  this  corresponds  to  a noise 
"power"  density  of  ,00204/( .1256 )^  = ,13  Y2/c/s.  Similar 
computations  may  be  performed  for  the  HP  filter.  The  equation 


corresponding  to  (5.5)  Is 


10.70  /V 


6.  Further  recommendations 

Consider  a simple  high-pass  filter 
0(a)  * s 

s + p 2 

operating  in  noise  having  a power  spectral  density  N /«  . 


(5.4) 


(5.5) 


(5.6) 


(6.1) 
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A simple  calculation  similar  to  equation  (f>.4)  shows  that  the 
r.m.s.  output  noise  level  from  G Is 

Thus  doubling  the  low-frequency  cutoff  from  p to  2p  reduces 
n by  a factor  of  /T  , or  3 db.  Because  the  low-  and  high- 
frequency  cutoffs  of  the  filters  discussed  In  this  beport 
(equations  (5.1))  differ  by  better  than  1.3  decades,  the 

A A 

variation  in  r.m.s.  output  noise  level  from  hlf  or  hhf  with 
moderate  variations  In  location  of  the  low-frequency  cutoff 
(.0213  or  .O89B)  will  be  approximately  as  given  in  equation 


(6.2). 


The  conjecture  was  made  in  section  I4.  that  frhe  low- 


and  high-frequency  characteristics  of  HLp  and  H^p  could  be 

designed  Independently.  While  this  is  not  strictly  true,  still 

It  Is  a useful  first  approximation.  Therefore  if  cmax  were 

halved,  from  I4I1..I4  to  22.2  sec.,  corresponding  to  a change  in 

the  the  original  specifications  of  seotion  1 to  8raax=750  ft.  or 

v . =UOkn.  (see  equation  (2.6)),  only  the  low-frequency  part 
ram 

of  HLp  (equation  ( 3 ' would  change;  the  principal  effect  on 

A 

Hlf  would  most  probably  be  simply  to  change  the  location  of  tbs 
lowest -frequency  pole  from  s*-.0213  to  ss-.Ol4.26.  As  noted  above, 
such  a ohange  would  result  in  a 3 db  deorease  in  r.m.s.  output 

A 

noise  level  from  HLp.  For  an  anomaly  signal  aotually  occurring 
at  the  high  end  of  ^p'8  bandpass  (say  at  t * 10.6),  the  change 
in  the  low-frequenoy  outoff  would  have  an  insignificant  effeot 
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on  the  peak  output  signal.  Hence  such  a signal  would  be  observed 
with  an  output  signal-to-noise  ratio  Increased  by  3 db. 

The  remarks  and  calculations  above  confirm  in  a quan- 
titative way  what  is  already  intuitively  known:  when  trying  to 
detect  a moderately  fast  but  possibly  weak  anomaly  signal,  an 
enormous  premium  is  paid  in  output  signal-to-nolse  ratio  for  the 
ability  to  see  slow,  but  very  strong,  signals. 

It  is  therefore  reoommended  that  the  low-frequenoy 
pole  of  HLp  be  adjustable  over  a range  in  absolute  value  from 
.0213  upwards  in  several  steps  to  perhaps  2 octaves  higher.  The 
adjustment  upward  from  .0213  would  be  made  by  the  MAD  equipment 
operator  whenever  it  is  known  that  his  target  is  moderately  near 
or  closing  moderately  fast.  There  seems  little  point  in  making 
adjustments  available  over  a wider  range  than  two  octaves,  be- 
cause H^p  would  then  be  approximately  designed  for  an  empty  set 

of  u values  (recall  Z r„  * Uii.U,  *c  , « 10.6,  and  two 

n&X|Lr  min* Lb 

octaves  from  Uii.U  is  11.1)*  That  is,  any  signal  observable  at 

A 

the  output  of  HLp  with  its  low-frequency  cutoff  more  than  two 
octaves  above  .0213  is  bound  also  to  be  seen  at  the  output  of 
ft 

"HP* 

It  is  possible  that  under  high  ambient  noise  conditions 
(large  values  of  Nc,  equation  (1.1))  the  output  from  HLP  would 

A 

become  totally  unusable,  and  even  the  output  of  HRp  would  be 
very  noisy.  Thus  it  seems  reasonable  to  provide  one  or  two 

A 

adjustable  low-frequency  eutoffs  for  Hnp  as  well;  it  is  expected. 
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however,  that  these  would  be  rarely  used. 

Swltchable  cutoffs  are,  of  course,  provided  on  present- 
day  airborne  MAD  equipment.  The  low-freouency  roll-off  of  these 
filters  is  much  faster  than  the  6 db/octave  of  the  filters 
proposed  in  this  report,  and  it  is  conjectured  that  the  reason 
for  this  is  to  minimize  noise  from  periodic  maneuvers  that  falls 
just  below  the  airborne  MAD  paasband.  For  shipborne  MAD,  periodio 

A A 

maneuver  noise  will  fall  within  the  passband  of  and  HHp, 
and  should  be  removed  by  notch  filters  as  recommended  in  section  1* 
It  is  not  recommended  that  the  high-frequenoy  cutoffs 

A A 

of  H^p  and  H^p  be  adjustable.  The  decrease  of  noise  power  den- 
sity with  increasing  frequency,  equation  (l.lV  Is  so  rapid  that 
the  resulting  increase  in  output  slgnal-to-noise  ratio  from 
such  provision  is  insignificant.  To  support  this  statement,  a 
filter  was  optimized  by  the  algorithm  in  Appendix  2 for  the 
single  signal  sQ(t;  UJ4..I4.).  It  would  be  expected  to  produce  a 
higher  output  signal-to-nolse  ratio  (with  Sg(t;  as  input) 

than  &Lp  with  its  high-frequency  outoff  moved  downwards  in  fre- 
quency. Yet  the  filter  optimized  for  s0(t;  yields  a 

signal-to-nolse  ratio  less  than  i db  better  than  unmodified. 
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Equations  (3*4)  were  derived  by  first  normalizing  the 
(bandpass)  spectra  of  aQ,  s^  and  Sg  so  that  the  maxima  of  the 
spectra  were  all  the  same,  and  Independent  of  c.  Then  at  each 
<■>,  Re {h ( j«e)}  was  taken  to  be  the  supremum  (over  aQ,  a^,  a 2 and 
t)  of  these  normalized  spectra.  Since  for  any  fixed  Xt  sQ 
occupies  a wider  frequency  band  than  either  s^  or  Sg,  equations 
(3.4)  result.  This  procedure  is  a sort  of  "worst  oase"  design 
that  ensures  H(J«o)  is  wide  efaough  to  pass  every  anomaly  signal 
under  consideration. 

Since  no  information  is  available  about  the  relative 
mean-square  values  of  the  coefficients  a0,  a^  and  ag,  it  18 
not  unreasonable  to  permit  the  design  to  be  dominated  by  the 
widest  band  signal,  a0. 

The  treatment  of  the  unknown  time-scale  o.d.f.  G(r) 
is  perhaps  less  satisfactory.  It  might  be  thought  that  without 
a priori  knowledge  of  0 that  the  so-oalled  "principle  of 
minimum  prejudice"  (PMP)  should  be  invoked.  The  PMP  asserts 
that  the  Minimum  prejudice  is  Inserted  into  the  problem  by  the 
designer  if  G(t)  is  chosen  to  be  the  distribution  with  the 
maximum  possible  entropy 

HCO]  * Jg(x)  log  g(r)  d*c 

consistent  with  the  constraints  of  the  problem,  where  g(r) 

■ dG ( r) /dr.  Por  the  problem  at  hand,  the  result  would  be  that 
Z should  be  assumed  uniformly  distributed  over  the  Interval 


(Al.l) 
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(2 .5>4,  Uli .U) . It  was  felt  that  this  ohoioe  for  G(r)  is  In- 
appropriate, and  the  PMP  is  not  applicable  so  simply.  Assuming 
■c  to  be  uniformly  distributed  would  certainly  resplt  in  a 
narrower  filter  than  the  filter  specified  by  equations  (3.4)# 
and  its  performance  on  signals  with e xtreme  values  of  c would 
necessarily  be  degraded.  Use  of  the  PMP  implies  that  some 
measure  of  filter  performance  Is  averaged  ober  all  values  of  c; 
poor  performance  at  extremes  of  the  range  is  more  than  compen- 
sated by  Improvements  at  moderate  r values.  But  the  design 
specifications  of  section  1 do  not  admit  such  an  averaging 
process;  the  filter  is  to  be  "good"  for  eaoh  value  of  r indiv- 
idually in  the  Interval  (2.54#  44.4) • 

The  choice  of  equations  (3.4)  to  specify  RefH]  is  by 
no  means  unique#  and  its  optimal  properties#  if  any#  are  ill- 
defined.  Nonetheless  it  has  the  virtue  of  simplicity  in  its 
formulation  and  in  the  final  filter  design  to  whloh  it  leads. 
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Appendix  4 - Glossary  of  Terms 


N(u>) 

s(t) 

s(t) 

s0(t) . Sl  (t),  Sa(t) 
30(t).  Si  (t),  Sa  (t) 
T 

S(ui) 

V(u>) 

H(juj) 


Kj  {• ) 

H(ju>) 


- noise  power  spectrum 

- anomaly  signal  of  a submarine 

- time  derivative  of  s(t) 

- three  classes  of  submarine  signals  comprising  s(t) 

- time  derivatives  of  sQ(t),  Si  (t),  sa(t) 

- ratio  of  smallest  slant  range  to  closing  speed 

- signal  power  spectrum 

- observed  signal  power  spectrum  (signal  plus  noise) 

- optimum  linear  filter 

- modified  Bessel  functions  of  the  second  kind 

- lumped  parameter  approximation  to  H(jiu) 
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cascade  with  the  basic  filters  in  order  to  remove  quasi-  cascade  with  the  basic  filters  in  order  to  remove  quasi- 
periodic  noise,  both  that  of  geomagnetic  origin  as  well  periodic  noise,  both  that  of  geomagnetic  origin  as  well 
as  periodic  maneuver  noise.  as  periodic  maneuver  noise. 
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From:  Commanding  Officer  and  Director 

To:  Commander,  Naval  Ship  Systems  Command  (SHIPS  1622F) 


Subj:  MEL  Report  537/66,  Transmittal  of 


1.  Transmitted  herewith  is  MEL  Report  537/66,  Linear  Filtering  for 
Shipboard  MAD  Recorders  (U),  NAVSHIPS  Sub-project  S-F101  03  21, 
Task  1522. 


2.  NAVSHIPS  Sub-project  S-F101  03  21,  Task  1522,  entitled  "Ship- 
board ASW  Magnetometer  Systems,  " is  an  exploratory  development 
effort  aimed  at  determining  the  feasibility  of  placing  magnetometers 
aboard  high  speed  surface  craft  for  ASW  search  and  classification 
purposes. 

3.  At  the  beginning  of  the  project,  it  was  realized  that  the  closing 
speeds  between  the  search  vessel  and  the  target  were  radically  dif- 
ferent from  those  encountered  with  aircraft  MAD  gear,  thus  neces- 
sitating redesign  of  the  magnetometer  filter  system.  Approximate 
filters  were  readily  derived,  but  a more  precise  filter  was  required 
for  optimum  signal  processing  in  a geomagnetic  noise  environment. 

To  accomplish  the  detailed  analysis,  Navy  Contract  61533-2551-65  was 
let  to  Dr.  Stephen  S.  Wolff  at  Johns  Hopkins  University,  a recognized 
expert  in  the  signal  processing  field. 
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4.  The  enclosed  report  is  the  output  of  this  contract.  It  is  mathemat- 
ical in  nature  based  upon  reasonable  assumptions  on  physical  parameters 
which  will  later  be  confirmed  or  modified  by  analysis  of  geomagnetic 
noise  data  obtained  during  the  summer  of  1966  at  Wallops  Island, 

Virginia.  The  main  .conclusion  of  the  report,  i.  e.  , that  a dual  output, 
two  band-pass  filter  is  required,  has  been  evaluated  at  MEL  and  will 

be  incorporated  into  future  systems  provided  that  the  maneuver  noise 
encountered  aboard  the  craft  carrying  the  sensor  does  not  radically 
affect  filter  design. 

5.  This  letter  is  unclassified  upon  removal  of  Report  53^/66. 
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